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Abstract 

We compare recent theoretical determinations of weak matrix elements of the electroweak 
penguin operators and Qg- We pay special attention to the renormalization scheme dependence 
of these determinations as well as to the influence of higher dimension operators in the different 
approaches. 



1 Introduction 



There has been recent progress in understanding the bosonization of some of the four-quark opera- 
tors which appear in the effective Lagrangian of the Standard Model describing K physics The 
bosonization of the electroweak penguin operators 

Qt = ^{sLl^dL) ^ eq{qRlpiqB) and Qs = -12 ^ eq{sLqR){qBdL) (1-1) 

q—u.d,s q—u,d.s 

in particular, where Cq denote quark charges in units of the electric charge and summation over quark 
colour indices within brackets is understood, has been obtained by two different analytic methods. 
One of the methods uses the framework of the 1 /N^ expansion [|j ; the other one combines dispersion 
relations with phenomenological input |Q, |5). Both methods are well rooted in the underlying QCD 
theory, and therefore they are competitive with more standard non-perturbative techniques based on 
lattice QCD simulations. In fact, the existence of lattice QCD estimates of matrix elements of the 
Q^ and Qs operators ^, ||] allows for a detailed comparison between these different approaches. 
It is precisely this comparison which is the main concern of this letter. We shall also discuss, within 
the particular case of matrix elements of the Qj operator, the role of "higher dimension operators" in 
weak matrix elements estimates which has been recently raised in refs. pi] . 

The operator Qy emerges at the Mw scale from considering the so-called electroweak penguin 
diagrams. In the presence of the strong interactions, the renormalization group evolution of Qj 
from the scale Mw down to a scale ji < rric mixes this operator with the four-quark density-density 
operator Qs: among others. The Qg operator plays an important role in the phenomenology of CP 
violation because, in the effective four-quark Lagrangian, it appears modulated by a Wilson coefficient 
Cg(/x^) which has a sizeable imaginary part induced by the Kobayashi-Maskawa phase in the flavour 
mixing matrix of the underlying Electroweak Model. 

It is well known that the bosonization of the Q-j and Qg operators leads to a term with 
no derivatives in the low energy effective chiral Lagrangian. This Lagrangian generates \^S\ = 1 
transitions among the pseudoscalar Goldstone fields of the spontaneously broken SU{2>)l x SU{3)r 
symmetry of the QCD Lagrangian with three massless flavours u, d, s. It has the following form: 

4A5|=i = -^^K<iK;^htr({/AfVtg^J+h.c., (1.2) 

where U denotes the 3x3 matrix field which collects the octet of pseudoscalar Goldstone fields, 
Qr = diag[2/3, — 1/3, — 1/3] is the right-handed charge matrix associated with the electromagnetic 

couplings of the light quarks and A^^"* is the effective left-handed flavour matrix ( A^^'^M = 6. 
= 1, 2, 3). Under chiral rotations {Vl, Vr): 



U^VrUVI, Qr^VrQrV^, Af^^FiAfVi, (1.3) 



and the trace on the r.h.s. of Eq. (1.2) is an invariant. Actually, this is the only possible invariant 
which can generate \AS\ — 1 transitions in the Standard Model to order 0(Gf a) in the electroweak 
coupling and to 0{p'^) in the chiral expansion. With the normalization chosen in Eq. ( |1.2| ), the 
coupling constant h is dimensionless. This constant plays a major role in the phenomenological 
analysis of if — > tttt amplitudes. It is one of the basic couplings of the low energy effective electroweak 
Lagrangian of the Standard Model that one would like to evaluate. 

The 1/Nc expansion in QCD jl^, Q offers a specific non-perturbative framework for discussing 



the dynamics which governs low energy constants like h in Eq. (1.2). It was recently shown in ref. P| 
that the contribution to the constant h from the Qy and Qg, four-quark operators can be calculated 
to first non-trivial order in the 1/Nc expansion. 

At the theoretical level, the comparison between different approaches can be made by evaluating 
the same matrix elements of the operators Qy and Qs in each approach; provided of course that the 

^See refs. |lj and Q and references therein. 
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calculations are made at the same scale ji and within the same renormalization scheme. In this respect 
it is important to remember, as clearly emphasized in ref. jist , that it is not enough to just state the 
definition of 75 in D dimensions, i.e. whether one uses 't Hooft-Veltman (HV) or naive dimensional 
regularization (NDR) or any other definition. One must also define the so-called evanescent operators 



and in this work we stick to the treatment and convention of ref. 1 16 . In particular we emphasize 
that within this convention Fierz symmetry is a valid concept even in D dimensions. 

However, there is still a technical difficulty when comparing theoretical evaluations of matrix 
elements. The calculations within the dispersive approach reported in ref. although aimed at 
evaluating matrix elements at 0{p^) in the chiral expansion, are not formulated within the large- 
Nc framework; while the lattice evaluations cannot be restricted, at least at present, to specific 
orders neither in the chiral expansion nor in the \/Nc expansion. Of course, if, as often claimed, the 
approximations within each approach are good, then the results should agree with each other within 
the estimated theoretical errors. Regretfully, the comparisons which have been made so far in the 
literature correspond to estimates which were obtained in different schemes. 

As seen within the combined frameworks of the chiral expansion and the 1/Nc expansion, the 
counting of possible contributions to matrix elements of the an d Q s, operators can be summarized 
as follows. Denoting by {Q7,i)\o(p°) the contribution to h in Eq. (|l.2|) from Qy.s, one has 



(Q7>|o(pO) =0(Ag+0(iV,"), and (Qg) b(pO) = 0(A^ + 0(iV,") , (1.4) 

were only the underlined contributions were calculated in ref. The contribution of 0{N^) to 
{Qfi)\o{p°) is Zweig suppressed. As we shall see, it involves the sector of scalar (pseudoscalar) Green's 
functions where it is hinted from various phenomenological sources that the restriction to just the 
leading large-A'c contribution may not always be a good approximation. It is precisely this issue 
which prevented us in ref. |^ from giving numerical values of the matrix elements of the Q% operator. 
Here we shall follow a different strategy that docs not suffer from the above shortcomings, which will 
allow us to make predictions for the matrix elements of the Qg operator as well. 



2 Bosonization of Qi and Qg 

2.1 The Q-j Operator 

Because of its left-right vector structure, the factorized component of the operator Q7, which is 0{N^), 
cannot contribute to order 0{p'^) in the low-energy effective Lagrangian. The first contribution from 
this operator of chiral 0{p°) is at least of 0{Nc) in the 1/Nc expansion and, formally, it can be written 
as follows 

(2.1) 
(2.2) 

(2.3) 



Q7 = 6 (Oi(^)) tr (u{x)Xf''^U^{x)Q 



where 

(Oi> = (0|(s-L7'^di)(Jfl7M5fl)|0) . 
Matrix elements of arc then given by the matrix elements of the effective operator 



tr ([/(x)Af )c/t(.)g^) = ^K-.+ + ^ 



on the r.h.s. of Eq. ( |2.l| ) times the dynamical factor 6(Oi). In particular, defining the isospin 
decomposition as (a = 7, 8) 



{TT+7r-\Qa\K^ 



I ((7r7r)/=o|Qa|i^° 



V2 



((7r7r)/^2|Qa|if°) 



(7r"7r"|QjX0) = I ((7r7r),=o|Qa|A'0) - z V2 {{^^)i=2\Qa\K'>) , 



(2.4) 



and taking into account that the operator of Eq. (2.3) has no transition between a K'^ and a two-TT*^ 
state, we find that 



2 



= -4(Oi), (2.6) 

where Q\ denotes the isospin 7 = 3/2 component of the Q^ operator, i.e. 

g(^/2' = 2sLl^dL {uB^lf^UR - d^f^dji) + 2sLj^'uLUR^f,dR . (2.7) 



Notice that the relations in Eq. (2^) are exact symmetry properties of the 0(p°) term in the chiral 
expansion. 

2.2 The Qs Operator 

As pointed out in ref. independently of large-iVc considerations, the bosonization of the Qs operator 
to 0{p^) in the chiral expansion can also be related to four-quark condensates by current algebra Ward 
identities, with the result 

Qs = -12 (02(a*)) tr (u{x)X^^^'>uHx)Qb)^ , (2.8) 

where 

(O2) = (0|(sLSfl,)(dfldL)|0). (2.9) 
From this result, recalling Eqs. (|2.4|), there follows the corresponding matrix element relations 



{{nn)j=,\Qs\K') = {{7:7:)j^,\Q'i^'^\K") = ^{n+\Q^^^'^\K+) = -l{n+\Qs\K+) (2.10) 
= §3(02), (2.11) 

where, again, the relations in Eq. ( p.lC| ) are exact symmetry properties of the 0{p^) term. The 

f3/2) 

operator Qg ' ' denotes the isospin 7 = 3/2 component of the Qs operator, i.e. 



1(3/2) 



= -4:{sLUR){uRdL) + i{sLdR) [{dndL) - {urul)] • (2.12) 



3 The large- A^c approach 

3.1 Evaluation of (Oi) 

As shown in ref. |^ the vacuum expectation value (Oi) can be expressed as an integral 

{0,) = l{-3^9,.J (04n-(,)) 
involving the two-point function 

nr^(g) = 2z / rf4^e^«-(0|T(L^(a;)7?''(0)t) |0) , (3.2) 



(3.1) 



with currents 

L^=q^{x)r^^qJ{^) and R^^ = q,(x)^'^l^ q^{x) , (3.3) 
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and i and j fixed flavour indices, i ^ j- This is the same two-point function which governs the 
electroweak tt"*" — 7r° mass difference |l7| ; however, because of the factorization of the short-distance 
contributions in the Wilson coefhcient Ct{h'^), the integral that appears in Eq. ( |3.lD is divergent for 
large . Then, for consistency of the matching between the long-distance evaluation which we are 
concerned with here, and the short-distance evaluation, the integral in Eq. (3.1) should be done in 
the same renormalization scheme as the calculation of the Wilson coefficients. The fact that we are 
interested in 0{p^) terms in the chiral Lagrangian implies furthermore that the calculation must be 
done in the chiral limit. 

As a first step, we shall therefore use dimensional regularization with Z? = 4 — e, and define the 
integral in Eq. (3.1) as follows 



where = —q^ and 



2\e/2 



D 



3(1 - D) (47r^^) 
167r2 r(2-e/2) 



dQHQ 



{~Q^Ulr{Q^)) \d , (3.4) 



The matching to short distances is controlled by the operator product expansion (OPE) ||T6 
D dimensions and in the large-iVj, limit one has 



hm (1 - D)Q^ {-Q'Ulr{Q')) 



D 



+ Oial)] ! + («:- 2/3)- (^V)' , 



where k depends on the renormalization scheme |^: 

K = -l/2 inNDR, and k = +3/2 inHV. 



(3.5) 
and in 

(3.6) 
(3.7) 



The next step is the use of the hadronic representation of the spectral function associated with 
^LniQ^), and this is where the large- A^c limit plays an important simplifying role. In this limit the 
spectral function consists of the difference between an infinite number of narrow vector states and an 
infinite number of narrow axial-vector states, together with the Goldstone pion pole: 



Mi 



(3.8) 



This spectral function is furthermore constrained by the fact that there are no operators of dimension 
d — 2 and d = 4 which contribute to the OPE of IlLRiQ^), which in turn implies the two Weinberg 
sum rules 

^fyM'y-Y.fAMl = F^ and ^ /^A^ - ^ /M = . (3.9) 

V A V A 

From the fact that Ii-LR{Q^) obeys an unsubtracted dispersion relation it then follows that 



- Q^^LRiQ^) = J2 



PaM\ 



QHQ' 



Ml) 



E 



f^M^ 



QHQ^ + M^) 



(3.10) 



Matching now the OPE of the function I1lr{Q^)\d in the quark-gluon language, as given in Eqs. (3.6) 
and ( [3.7[ ), onto the corresponding one in the hadronic language in Eq. ( 3.1Cl| ) requires 



Y^fyM^y-Y^flMl 



X + 



(3.11) 



D=4-e 



D=A 



The integral in Eq. (3.4) can now be made in the hadronic representation of large-Ac QCD. Using 
the MS-renormalization scheme, we get the following result 



MS "^^TT^ 



M2 



(3.12) 



^This is in agreement with the calculation reported in ref. 
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where 



= /i^ exp(l/3 



and K has been defined in Eq. (3.7). Notice that in ref. 
NDR and HV, in which k = 0. 



k), (3.13) 
01 another scheme was used, different from 



3.2 Evaluation of (O2) 

In the 1/Nc expansion one has that 

{02)^l[mnr + OiN'J] , (3.14) 

and therefore {O2) is directly related to the quark condensate {'tpip), provided we restrict ourselves to 
the first term in this expression which is 0{N^). However, as already mentioned in the introduction, 
there are reasons to suspect that subleading terms in the 1/Nc expansion involving some Green's 
functions with scalar (and pseudoscalar) density operators might be important, unlike those which only 
involve vector (and axial- vector) currents |^. Consequently one should be cautious about neglecting 
the 0{N^) contributions in Eq. ( |3.14 ) and consider instead 



(O2) = + {{sLSR){dRdL))c ■ (3.15) 

The unfactorized contribution involves Feynman diagrams which require gluon exchanges between 
at least two quark loops. These are the so called Zweig-suppressed contributions, which are indeed 
0{N^) in the 1/Nc expansion. This "subleading" term is governed by a two-point function of the 
type (Q^ = — , , and i, j = u,d,s with i 7^ j) 

^,,iQ')^z J dW'^-(0|r|q-(2;)i±^g,(x)g-(0)i^g,(0)||0). (3.16) 

More precisely 

{isLSH)idHdL))c = -( [ ^asm] , (3.17) 

where the integral should again be defined in the same renormalization scheme as the short-distance 
calculation of the Wilson coefficients; which explains the meaning of the MS subscript. There is, in 
particular, a contribution to this integral from the singlet 77'"' pseudoscalar, which in the chiral limit 
acquires a mass because of the axial U{1) anomaly. This contribution, which very likely plays an 



important role in the low-Q^ regime of the integral in Eq. (3.17), can easily be calculated with the 
result 



in agreement with results reported in ref. pl| , indicating the existence of 0{N^) terms which could 
be potentially important. In this respect we wish to point out that the failure to properly incorporate 
0{Nc) Zweig-suppressed contributions is also a problem which may seriously affect the quenched 
lattice calculations of matrix elements of Q% . 



To allow for the possibility of large deviations from the naive factorization in Eq. ( 3.14 ), it has 
become conventional to assume an ansatz of the type ||2^ 

(02)(/i)= ^ (V5V)'(m) , (3.19) 



^There are several phenomenological examples of this: then' mass, the possible existence of a broad a meson, large 
final state interactions in states with J = and / = 0, etc... [ |19[ In order to study this issue in a systematic way we are 
considering at present the possibility that the appropriate expansion for these exceptional Green's functions could be a 
1/Nc expansion in which rif/Nc is held fi xed where nf denotes the number of light flavours. This kind of expansion 
was originally advocated by G. Veneziano |20[. 
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where p{fj,) parametrizes possible deviations from the leading 0{N^) factorization where p ^ 1. 
The crucial observation B is that, on the one hand, the same vev (O2) also appears in the OPE 
which governs the high-Q^ behaviour of the I1lr{Q^) function (at = 4 and to lowest order in the 
perturbative evaluation of the Wilson coefficient |lq|). 



(3.20) 



On the other hand, we know that large- A^c QCD gives a reliable description for the IlLii{Q^) function 
since it involves the vector and axial- vector channels only. In fact, we can use the renormalization 
group to resum this lowest-order result into 0] 



lim I 



4^ 



p{Q) {4^4') 



-PiasiQ) 



01 



(3.21) 



where we have neglected the term proportional to (Oi), an approximation which we shall justify 

a posteriori. Here, {^ip) = {^'>P) ifJ-) i^og fi/ Ajjg)^^^^^ denotes the scale invariant quark condensate 
and p{Q) also includes the effect of next-to-leading perturbative corrections. Inserting the large-A'c 
expression which follows from Eq. ( 3.10| ) in the l.h.s. of Eq. (3.21) results (at large Q^) in 



A 



4^ 



p{Q) {^^) 



-PiasiQ) 



and therefore 



(02)(m) 



(3.22) 



(3.23) 



of Eq. ( |3^ ) 

> 1 as ^ oc 



due to the near 
. Since the l.h.s. 



It is worth noticing the remarkable independence of the r.h.s 
cancellation of the exponent = 2/11 and the fact that p{Q) 

is clearly Q independent, this is obviously a welcome result. 

We would like to pause here and comment on the 0{a^) term in Eq. ( 3.20 ). It is our understanding 
that the only existing calculation in the literature of this term [^3| was done in a different scheme 
of evanescent operators from the one of ref. Since this 0{a^) term is sensitive to the precise 

definition of evanescent operators ^ and since matrix elements of Qj^s eventually will have to be 
used in conjunction with Wilson coefficients to obtain a physical result, it is of the utmost importance 
that the extraction of (O2) in Eq. ( 3.20| ) be done within the same convention for evanescent operators 
as that employed in the calculation of Wilson coefficients. This is why we think, unlike the authors 



of refs. |j, that it is misleading to use the existing calculation of the 0{a^) term in Eq. ( 3.20| ) 



and, consequently, we refrain from including it to obtain (O2) in Eq. ( |3.23 ). This fact immediately 
introduces a major difference with respect to refs. || as the contribution of 0{a'^) they take turns 
out to be large, namely ~ 50%. 

Numerical estimates for (Oi) in Eq. (3.12) and (O2) in Eq. (3.23) will be made in the last section 
using the minimal hadronic ansatz approximation to large-iVc QCD (MHA). This approximation 
consists in limiting the large-iV^ spectrum of narrow states to the minimal number required to satisfy 
the OPE constraints relevant to the process that one is considering. In our case, this requires the 
consideration of one vector state and one axial-vector state, besides the pion pole 0. The MHA 
approximation has recently been shown p6|] to successfully reproduce sum rules which use the LEP 
experimental data of the ALEP H collaboration on hadronic r decay. This leads to a contribution of 
(Oi) to the r.h.s. of Eq. ( 3.2C ) which is indeed much smaller than the one obtained from (O2) in 
either scheme NDR or HV, justifying a posteriori the initial approximation which was made to obtain 
Eq. (H). 



^The analysis of ref. |24| shows that this term of 0{a^) can actually be reduced in half by a redefinition of the 
operator basis in the A — D extra dimensions. 

^ Sometimes we have referred to this particular case as the lowest meson dominance approximation (LMD) tea]. 
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4 The dispersive approach and the role of "higher dimension operators" 



When computing (Oi) through Eq. (3.1), the authors of ref. |4| have chosen to spht the integral in the 
regions < < fJ,^ and fj,^ < < oo where fi plays the role of a sharp momentum cutoff, i.e. not to 
be confused with the scale appearing in dimensional regularization. The low-Q^ part is amenable to 
the use of the experimental data through a dispersion relation. The high-Q^ part, on the other hand, 
is divergent and must be renormalized, e.g., in dimensional regularization in Z) = 4 — e dimensions 
with minimal subtraction (MS). In this way one obtains 

(Or)^(M) = (Oi)--«(m) + ^ 32,f tT^I^ /J ^^'(Q')'"^^' {^Q'llLniQ')) , (4-1) 
where (Oi)'^"*°'^(/i) corresponds to the low-Q^ part, i.e. 

\cutoff/,,N _ 3 /"°° ,2 ^ + 1 



(Oi)™*°**(a^) = / dtt'log—^-lmULRit) . (4.2) 

3z7r^ Jq t tt 

Notice that, when writing Eq. (^^), one is using the first and second Weinberg sum rules. 

We are now in the position to discuss the issue of "higher dimension operators" in the calculation 
of weak matrix elements which has been recently raised in refs. 0. On the one hand, these 
authors have stressed the fact that the long-distance evaluation of weak matrix elements should be 
done in the same renormalization scheme as the short-distance evaluation of the Wilson coefficients. 
Although this is a well known fact, they observe, quite rightly, that most of the available long-distance 
calculations of weak matrix elements are done in a cut-off scheme while using the MS results of the 
Wilson coefficients. Exceptions to that are the lattice QCD simulations and the large-TVc analytic 
results in refs. [Q, |^. On the other hand, a more subtle issue also raised in refs. is the fact 
that in doing integrals of Green's functions which govern the long-distance evaluation of weak matrix 
elements, in the same MS renormalization as used for the Wilson coefScients, one is still confronted 
with the fact that, in principle, the full string of terms in the OPE contribute to these integrals because 
the integration over the euclidean momentum goes all the way to infinity. We entirely agree with this 
observation but we want to show how the large-iVc approach we have been advocating avoids this 
criticism while the dispersive approach, as used in ref. Q|, fails to incorporate the effect of higher 
dimension operators. 

We consider first the evaluation of the vev (Oi)*^'^ in Eq. (4.1) to illustrate the point. Depending 



on the input for Ulr in the integral over the high-Q^ region in this equation, one may obtain a 
unsatisfactory determination of (Oi)^^ even if (Oi)'^"*°^ is derived reliably from the experimental 
data. The authors of ref. M input for Ulr its leading asymptotic behaviour coming from dimension 



six operators in the OPE, i.e. the equivalent to our large- A^c Eq. (3.6). While for /i large enough this 
is certainly a reasonable thing to do, for the realistic value oi fi — 2 GeV the question still remains 
as to how large the contribution from higher dimension operators in the OPE may actually get to 
be 0,|ll| 



hi this regard we would like to point out that , in the large-iV^ approach, the integrand in Eq. ( 3 . 1 C ) , 
if expanded in powers of gives not only the leading power of the OPE, which we use as a 

constraint, but also the higher powers. As a matter of fact, and this is the crucial point, in ref. |Q 
it is shown that the MHA approximation to large-A"c discussed at the end of section (||) does a 
remarkably good job in predicting these higher vev's as compared to the ALEPH data. Therefore 
one can use the MHA approximation to confidently evaluate the size of the contribution from higher 
dimension operators to the high-Q^ region in the integral of Eq. (4.1). One then obtains 



(Oi)^^^(m) - (Oi)^"'°'^(m) - |^«s(02) [l + n] (4.3) 



' ^E/v^^vlog^^^^-E/MloST^l ' (4-4) 
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where the expression (4.4) reads in the MHA approximation {gA — My/M\) 



327r2 1 ~ 



log- 



— log 

9A 



9A 



^1- 



(4.5) 



Equation (O) is the result of ref. ]4| (in our notation), with only the contribution from dimension-six 
operators. The expression in Eq. ( [4.5| ), which represents the contribution from operators of dimension 
8 and higher, indeed goes to zero as an inverse power for large values of fi, as expected. However 
it amounts to a ~ 50% reduction at /i = 2 GeV in NDR (see next section). We emphasize that in 
ref. 1^ the MHA approximation was shown to yield correct predictions in the OPE up to operators 
of dimension 10, so we consider Eq. (4.5) to be reliable at least up t o thi s order. 

The corresponding discussion in the case of (O2) starts at Eq. (3.20) with the contribution from 
(Oi) neglected, as already explained. Therefore (O2) is the operator modulating the l/Q^ fall-off in 
the OPE of the function Ulr- At large Nc one obtains 



(O2 



(4.6) 



where 



1 - 



1 - 



Ml 



MP 



Ml 



9a(. 



Ml 



(4.7) 



where the last expression is a consequence of the MHA approximation. Again C{Q^) — 1 vanishes as 
an inverse power at high Q^. Since in refs. ^, p| the contribution from higher dimension operators 
is neglected, these authors effectively take C{Q^ = 1. However a typical value for C{Q^) is ^ 40% 
above unity for Q = 2 GeV (see next section for details). 



5 Numerical Estimates and Comparisons 

Using Eqs. ( p.l2[ ), ( p.l3| ), ( |3^ ) and ( p.23[ ) one obtains in the MHA approximation to large- A^^c QCD 
the simple expressions 



(01) (A*) = 

(02) (m) 



3 F^Ml 



32^2 i^g^ 
1 



log 



9A 



(9_A^ 
\Ml 



FiMl 



16nasi^l) gA 

For their numerical evaluation wc shall take ]2^ 

My = (750 ± 25) MeV , i^o = (87 ± 3) MeV , 5^ = 0.50 ± 0.06 , 
and as{fJ. — 2GeV) — 0.33 ± 0.04 . Then the previous expressions lead to 



(01) (^=2Gcv) 

(02) (^=2Gcv) 



(-1.9 ± 0.2) X 10-5 GeV*^ in NDR 
(-1.1 ± 0.2) X 10-4 GeV*^ in HV 

(2.9 ± 0.6) X lO-'^ GeV^ both in NDR and HV 



(5.1) 
(5.2) 

(5.3) 

(5.4) 
(5.5) 
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which in turn imply, through Eqs. (|U), ( px| ) and 

in l;<'0\ (0.11 ± 0.01) GeV^ in NDR , 

((^^)/=2|Q7|A )(^=2G.v)= (0.67 ± 0.09) GeV^ in HV ^ ' ^^'^^ 

and 



((^^)/=2|g8|if")(,,=2Gev) = (3.5 ± 0.8) GeV^ . (5.7) 

In order to facihtate the direct comparison with the lattice results we also give the equivalent 
result for the matrix elements 



/ +in(3/2)|i^+\ _ (-9.8 ± 0.6) X 10-3 GeV^ in NDR , 

\Qr \K )(^=2ocv)- („5.8±o.8)xlO-2GcV4 in HV ^^■^> 



and 

(^+|gf/''|i^+)(^=2Gcv) = (-0.30 ± 0.07) GeV^ . (5.9) 

These are our predictions for the O1.2 condensates and for the Qt^s matrix elements. The errors quoted 
come only from the propagation from the input values for Afy, Fq, qa and as(2GeV). We remark 
that, unlike the case of Q7, matrix elements of Qs do not show any dependence on the renormalization 
scheme (NDR vs. HV) at this level. 

In Table 1 we give a joint comparison of our results with the existing different evaluations of matrix 
elements with which we can compare scheme dependences explicitly ^. 



Table 1: Summary of matrix elements My^g = ((7r7r)/=2|Q7.8|^'^)(2Gev) using naive dimensional regu- 
larization (NDR) and the 't Hooft-Veltman scheme (HV), in units of GeV'^. 



Matrix Elements 


M7(NDR) Af7(HV) Af8(NDR) Af8(HV) 


rcfs. [' 


i> 'A 


0.22 ±0.05 1.3 ±0.3 


rcf. |ei 


0.11 ±0.04 0.18 ±0.06 0.51 ±0.10 0.62 ±0.12 


This work (see also ref. Q) 


0.11± 0.03 0.67±0.20 3.5±1.1 3.5±1.1 



In the table we have rounded off the errors of our predictions to an overall 30%, which we believe to 
be a generous estimate of the systematic errors in our approach |3(]||. Fo r the case of Mg this has the 
caveat that it of course assumes that the 0{a1) corrections to Eq. ( ^.2()| ), once computed in the right 
operator basis jl^, will be of a reasonable size. 

We find that, within the combined errors, our results for Mi are in agreement with the lattice 
results in the NDR scheme, but not in the HV scheme while they disagree by a factor of two with 
the dispersive results. The origin of this disagreement can in fact be traced back to the effect of higher 
terms in the OPE found in Eq. (4.5) which, as discussed in the previous section, have been ignored 
in the dispersive approach. As to the results of Mg, we are in disagreement with both the lattice 
and the dispersive results. The discrepancy with the former may originate in the fact that most of 
the contribution comes from an OZI-violating Green's function which is something inaccessible in the 
quenched approximation. As a matter of fact, lattice results are compatible with a value p ~ 1 in 
Eq. (3.19) whereas we are finding that p ~ 6 for (7/j-0) (2 GeV) ~ (—0.240 GeV)^, which is also the 
phenomenological result found in rcf. On the other hand, we do agree with the lattice results on 
the independence of scheme (NDR vs. HV) in Mg. As to the discrepancy on Mg with the dispersive 
results, it can be traced back to the ~ 50% 0{al) correction they use for the Wilson coefficient in 

^This explains why, in this table, we only quote the lattice results of ref. Model dependent calculations of 

the so called B factm-s associated with the Qy^g operators can also be found, in order of increasing sophistication, in 



refs. |g8|, |2l| and 

'^Although there is a difference bewteen what the Rome and Munich groups call the "HV" scheme lla], we have 
checked that the values for My(HV) quoted in the table are indeed consistent with the ones obtained fromMT {N D R) 
and Ms (NDR) and the mixing of Q7 and Qg under the appropriate change of scheme; therefore, the disagreement with 
the lattice result for My[HV) is correlated to the strong discrepancy we have with the lattice result for Ms{NDR). 
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Eq. (Og) plus the ~ 40% effect coming from the contribution from higher dimension operators in 
Eq. 
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